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The present invention relates to methods and 
apparatus for analysing quasi-repetitive signals 
con5>rising a series of similar complexes* and in 
particular, but not exclusively, to such methods and 
apparatus suitable for analysing heart beat signals such 
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The condition and behavio\ir of the heart can be 
measured in a variety of ways. Electrical behaviour can 
be measured using electrocardiography (BCG) . Electrodes 
are placed on a subject's abdomen and small electrical 
signals whidh relate to corresponding electrical signals 
in the heart are recorded. Electrodes, especially 
reference electrodes, may also be placed elsevdiere on the 
body. Placement of multiple electrodes at selected points 
on the body enables a variety of different signals to be 
picked up simultaneously (multiple lead 
electrocardiogra£»hy) , from which a reasonably complete 
picture of the electrical activity of the heart, and any 
related clinical abnoimalities can be deduced. A heart 
beat con^jlex from a typical EC<3 signal is shown in 
25 sin5?lified form in figure 1. The signal coinprises many 
such coiiplexes in series, each corresponding to one heart 
beat. An ECG signal may be taken from each lead. 
Typically^ tts^elve leads ere used. 

30 The principal features of an ecg con©lex have widely 

accepted designations. An initial P-wave is followed by 
a larger QRS-wave. A T-wave follows the QRS-wave. Various 
minor variations in the shape of a complex, as well as 
missed and extra beats are ea^ected in a healthy subject. 

35 Heart abnormalities, diseases, and drugs also affect the 
shape and timing of the feattures of a complex. Trials of 
a drug generally include experiments to determine any such 
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changes, from which the effect of a drug on the 
electrocardiographic activity of the heart can be deduced. 

Pre-clinical drug trials generally include 
administering a drug to animals using various dosage 
sc^hemes. During an experiment, tarhich typically lasts 
several hours, one or several doses of a drug may be 
administered and an ECO recording made. It is time 
consuming to manually analyse the large amount of ECG data 
resulting from such escperiments , and from other situations 
in which long series of heart beats are recorded. 



Automatic analysis of BCG complexes to identify 
various parameters is well known. Analysis of the mutual 

15 timings of various features of signals from one or more 
ECG leads, for exaitrple, is routinely built into EOG 
machines for clinical use. Such shape analysis is usually 
carried out on single selected heart beat consplexes. 
Automatic shape analysis of the average of a series of 

20 complexes is also known, for example from OS 5# 139^027, 
which describes an apparatus which also automatically 
looks for features in an BCG signal such as narrow spikes, 
electrical noise and other artifacts and discards them 
before forming an average of a plurality of complexes. 

2S The averaging of a plurality of complexes before 

performing shape analysis is advantageous in that it 
reduces the computational cost of the expensive shape 
analysis process and improves the signal to noise ratio of 

30 

However, the known schemes for automatic analysis of 
BCG signals are not ideal in many situations, such as pre- 
clinical drug trials, because they do not adeqpxately allow 
for natural variations in shape and timing between the 
35 members of a group of heart beat complexes to be averaged. 

The invention seeks to address the above and other 
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problems of thja related prior art. 

Accordingly, the inv«Jtion provides a method of 
analysing a heart beat signal which can^rises a series of 
heart beat complexes, the method conpriBing the stetps of: 

selecting a set of said complexes; 

identifying from the set a subset of coinplexes each 
of which satisfy a similarity criterion with respect to 
members of the set; 

calculating a representative conplex from the ta 
of the subset; and 

establishing one or more heart beat signal parameters 
from said r^resentative conplex. 



15 aaie heart beat signal is preferably an ecG signal. 

but the method may also be applied to other types of heart 
beat signal, such as a signal from a pulse oximeter, a 
thoracic motion sensor or an ultrasound device. Each 
heart beat coinple^ preferably corresponds to a sixigle 

20 heart beat, or portion thereof. Each selected set of 
eomplexee preferably includes only genuine conplexes. 
excluding artifacts of various sorts familiar to the 
Skilled person, such as spikes and electrical noise. The 
complexes in any one set are preferably consecutive within 

25 said signal. Typically, the set will coirqprise between 
about two and twenty complexes, and preferably between 
four and tan complexes. 



30 The heart beat signal parameters may be calculated 

from the representative conQ>lex using a variety of 
tedmiques known in the art. For an BC6 signal, 
parameters such as the PQ, qrs and QT intervals, and 
parameters derived from intervals between features on 

35 different leads may be calculated. 



01:00:91 eoo3*60*93 eio 

^la £08 'ON 



- 4 - 

Preferably^ the step of identifying comprises 
identifying from tlie set a subset qf coirplexes each of 
which satisfy a similarity criterion with respect to an 
average of the members of the set, where the average could 
5 be a mean, an lumormalised stun, or a variety of other 
types of average. Preferably, said step of calculating 
coxqprises calculating the average complex of the tnembers 
o£ the sxjibset, where again the average could be one of 
several types of average. 

10 

Preferably, the inethod includes a step of analysing 
said heart beat signal to identify and select genuine, or 
valid heart beat complexes to form the above mentioned 
series I before any of steps mentioned above are carried 
15 out* The method may also include steps of aligning the 
complexes of the series, the set and/or the subset prior 
to operations such as correlation and averaging 

preferably, the similarity criterion conqprises a 
20 threshold calculated using the congplexes of said set and 
an average of these complexes. For exajtrple, the criterion 
could be based on coefficients of correlation between each 
member of the set and the average of the set, the 
criterion comprising a threshold correlation coefficient. 
25 Such a threshold coefficient may itself be based on the 
distribution of the calculated correlation coefficients 
for the set. However, the similarity criteria preferably 
also coir^rises a predefined threshold value, for example a 
fixed correlation coefficient of 0.98, above which the 
30 criterion is always satisfied. This ensures that no 

complexes fail the similarity test if they are all very 
sixnilar. 

Said step of identifying may further comprise 
35 aligning the complexes of said set, forming a set- 
average complex from said aligned coxtvplexes, comparing 
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each congplex of said set to said set-average con^)!©*, 

deterroining if each result of comparing satisfies 
said similarity criterion, and forming said subset from 
only those coinplexes of said set for which the similarity 
criterion is satisfied. Preferably, the step of comparing 
conprises correlating each coniplex of said set with said 
set-average eoaj^lex. 

If said signal is an electrocardiograx>h signal, said 
step of aligning the complexes of said set preferably 
comprises the step of aligning said coxrplexes on an R-tpp 
feature of each coniplex, although other alignaient schemes 
could also be used. 



Vha method steps described above may conveniently be 
implemented using ccanputer software for execution on a 
suitable conputer system. Such software may conveniently 
be written on one or more cona?uter readable media such as 
CDEOis. The invention also provides a computcsr system 
arranged to carry out the methods. Such a computer system 
may provide, for example, an extract set element adapted 
to select a set of complexes from the series of complexes 
to be analysed, a select subset elem^t adapted to 
identify from the set a subset of complexes each of which 
25 satis^ a similarity criterion with respect to members of 
the set, and a coirtbiner element adapted to calculate a 
representative compless from the members of the subset. 

Otoe invention also provides apparatus for analysing a 
heart beat signal con^prising a plurality of heart beat 
complexes, the apparatus congprising: an analysis engine 
adapted to automatically calculate a representative heart 
beat complex for each of a plurality of preselected 
intervals of said signal; and 

an editor adapted to enable a user to edit calculation 
and/or heart beat parameters relating to each 
automatically calculated representative heart beat 
congplex. Typically, each interval contains a set of heart 
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beat coxnplexes. The analysis engine is preferably adapted 
to automatically decide which con5>lexes in the set to 
include in each representative heart beat coinplex^ and the 
editor is preferably adapted to enable the user to 
5 subsequently change which complexes are included. 

* 

' The editor niay enable the user to control calculation 

j parameters which detemdne whether a particular coi^plex 

should be included, not included, or included subject to 
10 an automatic or predefined similarity criterion, in the 
corresponding representative heart beat complex. Thie 
r€5)resentative complex may be an average of the included 
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The editor may further provide the user with a 
graphical display of at least some of the heart beat 
con^lexes within a current one of the intervals of saj 
signal, and with a graphical display of the current 
representative heart beat complex for the current 
interval. In this way, the user may be provided with 
facilities to validate the representative heart beat 
coniplex and correct it if necessaaty. 

For exaxnple, the editor may be further adapted to 
25 display, on the display of the cxxrrent representative 

heart beat complex, one or more feature markers indicating 
the locations of one or more automatically determined 
features of said representative heart beat conqplex, such 
as P, Q, R, S and T features of an ECG complex. The 
30 editor is preferably further adapted to enable the user to 
move said markers, for example by means of a drag and drop 
action using a pointer device. 

Preferably/ the analysis engine is adapted to 
35 calculate one or more heart beat parameters, especially 
shape parameters, of each representative heart beat 
complex, and said editor is adapted to display said heart 
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beat parameters and to update said display according to 
changes made by said user. Such heart beat parameters 
could include, for ecg analysis, the duration of the qt, 
PQ and QRS as %rall as of other intervals between features 
for a single lead or between leads. 

m 

The apparatus is preferably provided by means of a 
coinputer program and appropriate data elements, which may 
be provided on one or more removable eon^uter readable 
media, executing on a coniniter. One preferred embodiment 
is implemented using Microsoft Excel (RTM) spreadsheet 
templates, macros and other appropriate data structures 
and files. 



The method and apparatus may be used in a variety of 
circumstances and applications. In the embodiments set 
out below, pre-clinical drug trials are particularly 
mentioned, but other applications include any analysis of 
heart beat data including clinical drug trials, medical 
check UPS, patient monitoring in intensive care, during 
and after surgical procedures, in athletic activity 
monitors and so on. 

The methods and apparatus may also be applied to 
repetitive or opiasl -repetitive signals other than heart 
25 beat signals including to signals of medical origin and 
use and to non-medical signals. 



ESribodiments of the invention will now be described, 
• by way of isxample only, with reference to the accompanying 
30 drawings, of which: 

Figure 1 shows features of a typical ECG signal; 

Figure 2 sets out a scheme for analysing a heart beat 
signal ; 

Figure 3 illustrates how steps 14 and 16 of figure 2 
35 may be put into effect, in more detail; 

Figure d illustrates apparatus and/or coxuputer 
program elements for putting the method steps of fig\ires 2 
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and 3 into effect; 

Figure 5 sets out a scheme for analysis and 
validation of a heart beat signal; and 

Figure 6 shows a graphical user interface for 
5 assisting in the validation step of figure 5. 

Figure 2 illustrates, schematically, a method of 
analysing a heart beat signal 10, such as an ECG signal, 
eRibodying the invention. The method conqprises steps of 

10 selecting a set of heart beat complexes from within the 
signal (step 12) p identifying a subset of those coniplexes 
which all match a quality, similarity or mutual similarity 
criterion (step 14) and calculating an average or 
representative conqplex of the identified subset (step IS) - 

15 The average complex may be stored and used to calculate 
one or more heart beat parameters (step 18) . If the 
signal is an ECG signal, then these heart beat parameters 
may be, for exan«>le, a QRS interval, a PQ interval, a QT 
interval and so on- 

20 

St^s 14 and 16 of figure 2 may.be implemented as 
illustrated in figure 3. The members of the selected set 
of heart beat complexes 20 are first aligned in st^ 22 # 
The alignment may conveniently be carried out using a 

25 single dominant feature such as the R-peaJc of an fiCO 

cox(^lex, or may involve optimization of alignment of more 
then one feature or region of each conplex. After 
alignment, the con^lexes are then summed or averaged in 
step 24. A correlation coefficient between each aligned 

30 con^plex ajxd the averaged conqplex is then calculated in 

step 26. *Ehe set of correlation coefficients so obtained 
is used to calculate a threshold parameter in step 28. 

The threshold parameter calculated in step 28 is used 
35 in step 30 to discard one or more of the original 

complexes 20 according to a similarity criterion using the 
parameter. The remaining subset of coniplexes may again be 
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for exainple with reference to an R-top, in step 
32, and the subset is then used to calculate an average 
complex, representative of the original coinplexes 20, in 
step 34. 

A suitable threshold parameter calculated in step 28 
is the value above which 95% of the correlation 
coefficients would be expected to lie, based on the 
average and standard deviation of their actual values. 
The similarity criterion used in step 30 could then be 
such as to discard coitplexes having a correlation 
coefficient calculated in step 26 falling below this 
threshold parameter, subject to a maxiinum threshold 
parameter value of 0.98. 

The method illustrated in figtireg 2 and 3 inay be put 
into effect on any suitable conputer equipment, in 
particular using suitable software elements which may be 
provided on on© or more conK>uter readable media- Suitable 
elements of such software, or con^uter apparatus suitably 
programmed, are illustrated in figure 4, The arrangement 
may be controlled and the results of the iiftpl^ented 
process displayed using the apparatus and user dlnterface 
further described below in relation to figxire 6. 



Figure 4 illustrates a memory device 40, such as a 
disk drive or HAM, storing an ECG or other heart beat 
signal. The signal may be preprocessed in various ways by 
preprocess element 42- The preprocess element 42 may 
carry out cleaning and filtering processes on the signal, 
and in particular may identify and label valid heart beat 
complexes. Noise artifacts and badly distorted coinplexes 
are rejected, for exaj^ple by comparison against texuplates 
representing acceptable complex shapes, hy checking the 
timing of each conplex within the signal, or by other 
techniques known in the prior art. 
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An extract element 44 selects one or more sets o£ 
^ coirgplexe:s for further analysis t each set preferably 

I containing a fixed niuriber of consecutive complexes. For 

the purposes of pre-clinical drug trials each set could be 
5 a series of ten coiqplexee selected at fifteen minute 

intervals throughout an eacperiment , with extra set© being 
selected in each period iinraediately following a drug 
administration • 



10 An align element 46 aligns the complexes of each set* 

A calculate criterion elaiient 48 calculates the parameter 
or parameters, such as a threshold, for deciding which 
con^lexes of the set to use further, and which to reject » 
for example as described above in connection inrLth steps 

15 24, 26 and 28 of figure 3. An identify siabset element 50 
uses the parameter or parameters from the calculate 
criterion element 48 to select a subset of the con^lexes« 
as described above in coxonection with step 30 of figtire 3. 
Finally* the selected subset is used by averaging element 

20 52 to calculate a representative conplex* and the analyse 
conplex element 54 automatically analyses the average 
conplex to obtain desired heart beat parameters. At each 
stage of the process / and especially at the final analysis 
stage^ results may be written back to the memory 40 in a 

25 manner so to be suitably associated with the original 
signal data. 



Generally, and particularly for pre-clinical drug 
trials, the heart beat signal analysis scheme described 

30 above is incorporated into a process such as that set out 
in figure 5. An ECG signal is recorded at step 60. 
Corresponding experiment information is input at step 62 » 
and used to perforin an automatic analysis of the ECG 
signal at step 64. The results of the automatic analysis 

35 are checked, and if necessary adjusted at validation step 
66, and finally an experiment report is produced at step 
68. Ihe e^cperiment information may include conditions of 
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the escperimenh, such as time, date, duration, location, 
subject animal, drug, dosage and so on, as well as 
parameters expressly constraining the subsequent 
calculation step 64 such as the time points at ^ich 
detailed analysis is required, how that analysis is to be 
carried out, and what result heart beat parameters are to 



In one embodiment of the invention, steps 62 to 68 of 
fixture S are in5>lemented using Microsoft Excel (RTM) 
spreadsheet teinplates and macros, with ECG sigzial data 
being inserted in a conrpatible file format. Figure 6 
illustrates a graphical user interface provided as part of 
this embodiment to implement the validation and adjustment 
step 66 of figure S following the calculation step 64. 
'Phis interface is displayed ty a Kicrosof t KKcel iRTtH) 
tool executing on a conventional personal computer, 
providing conventional pointing device, keyboard and 
visual dis5play peripherals. 



The diiration of the ^cperiment being validated in the 
interface of figure 5 is shown in the time point area 70. 

The time points (in minutes) at wliich detailed analysis 
of the ECG signal is required are shown • The first time 
25 point is minus fifteen minutes, because the zero minute 
point is the time of a first drug administration. A 
current time point is indicated by a box 72 surrounding 
the ""IS" value. Completed \''alidatian at a time point is 
indicated by colour coding in the time point area 70- 



The ECG parameters being validated are based upon up 
to twelve time-parallel signals from up to twelve ECG 
leads, listed in the left hand coluimi of a validation 
table 76. A current ECG lead is indicated by a box 78 
35 surrounding the lead II row. Colour coding in the 

validation table 76 illustrates which beats of which lead 
have been validated, and which leads do not require 
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validation. A protocol editor may be available. In said 
protocol editor the user can predefine time points, leads # 
parameters and so on tbat are to toe used or calculated in 
the BCG editor. The horizontal axis of the lead table is 
5 representative of the ten beats (or other nuniber as 
appropriate) to be averaged at the current time point. 

The ECO signal of the beats to be used at the current 
time point for the current lead are plotted in a beat plot 

10 area 80, as multiple curves overlaid according to an 

automatic alignment process. The ten beats shown in the 
beat plot area 80 of figure 6 have been aligned on their 
R-top features. Which beats are plotted can be controlled 
using a beat plot selector 82, which currently indicates 

15 that all beats of the set except beats which have been 

rejected on the basis of a similarity test (ie used beats) 
are shown* Other options provided by selector 82 are all 
beats of the set, and single selected beats. 

20 Below the beat plot area is an average beat plot area 

86 which displays, using a horlrcontal axis aligned with 
that of the beat plot area 80, an average beat calculated 
according to the method illustrated in figures 2 to 4. 
Also shown are markers B8 (P.Q/S and T) showing the 

25 positions of corresponding features of the average beat 
which have been automatically calculated, for example 
using beat shape analysis techniques known in the prior 
art. The markers 88, which could mark various different 
features, can be moved by the user. Heart beat parameters 

30 derived from the current marker positions are given in a 
heart beat parameter area 90. 

Although the average beat and consequent heart beat 
parameters are automatically calculated for all time 
3S points according to a schedule included with the 

information provided in step 62 of figure 5, the user of 
the interface of figure 6 can subsequently validate and 
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: Intervene in the beat averaging process by means of the 

I controls provided in the average beat control area 94. 

! Each row of this area corresponds to one of the beats at 

the current time point. The "Correl' column shows the 
5 correlation coefficient of each beat, once aligned, with 

the average of all beats of the set, and the "Correl tJ- 

colxann shows the correlation of each beat, once aligned, 
i with the average of all beats currently used to form the 

average beat shown in the average beat plot area 8€ 

10 

The "User- column of the average beat control area 94 
enables the user to allow the tool to decide whether to 
include a beat in the average according to the 
predetermined scheme, or to override the scheme and choose 
to include or eatdude the beat from the average. iSiis 
facility allows a user to reject a distorted or otherwise 
incorrect or undesirable beat, and to include an otherwise 
rejected beat if necessary. 

Navigation area 98 provides facilities for moving the 
current validation point between leads, beats, time points 
and experiments, for accepting and validating data and 
derived parameters and so on. 

Alternatives and variations to the ertbodiments 
described above will be apparent to the person skilled in 
the art, while still falling within the scope of the 
appended claims and their equivalents. For example, a 
variety of different programming languages could be used 
to implement the invention, on a variety of different 
hardware platforms including dedicated BCO machines, 
portable heart beat monitors, and conventional computer 
work stations. A validation interface could include 
various combinations of graphical elements selected from 
those illustrated in figure 6, as well as other elements. 
The validation interface may also be provided as a 
separate entity to any automatic processing means, the 
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results from which it is to be used to validate. 

Varioud criteria may be used to select which heart 
beat complexes to include when calculating a 
5 representative complex:. In some applications only 

complexes deviating to a large extent from the average of 
a set will be rejected^ while in other applications only 
near identical coxr^lexes will be included. 
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1. A method of analysing a heart beat signal which 
5 conqprises a series of heart Jseat con5>le55eB, che 

method comprising the steps of : 
selecting a set of said complexes; 

identifying from che set a subset of con©lexes each 
of which satisfy a similarity criterion with respect to 
10 members of the set; 

calculating a representative coaBplex from the meaiibers 
of the subset; and 

establishing one or more heart beat signal parameters 
from said representative conqpleix. 



3S 



2. The method of claim 1 wherein said step of 
identifying comprises identifying from the set a subset of 
eomplexes each of which satisfy a similarity criterion 
with respect to an average of the members of the set. 

3. Ml© method of claim 2 wherein said similarity 
criterion comprises a threshold calculated using the 
complexes of said set and an average of the members of the 



4. The method of claim 3 wherein said similarity 
criterion further co«©rises a predefined threshold value 

5. The method of any preceding claim wherein said step 
of calculating a representative complex is a step of 
calculating an average complex of the members of the 



6. The method of any preceding claim further comprising 
a step of analysing said heart beat signal to identify and 
select genuine heart beat complexes to foxro said series. 
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7 . The method of any preceding claim wherein said step 
of identifying coTRprises: 

aligning the con^lexes of said set; 

forining a set-average complex from said aligned 
coxiiplexes ; 

coxttparing each complex of said eet to said set- 
average coniplex; 

determining if each result of conparing satisfies 

said similarity criteriOTi; and 

forming said subset from only those coitgplexes of said 
set for which the similarity criterion is satisfied. 



8. The method of claim 7 wherein said st^ of coxrparing 
con«>rises correlating each eoirplex of said set with said 

IS set-average complex^ and comparing each resulting 
correlation coefficient with a threshold value. 

9. The method of claim 7 or 8 wherein said signal is an 
electrocardiogram signal and wherein said step of aligning 
the complexes of said set conprises the step of aligning 
said conqplexes on an R-top feature of each couc^lex- 



10. The method of any of claims 1 to 8 wherein said heart 
beat signal is an electrocardiogram signals 

11. Coinputer apparatus adapted to carry out the method 
steps of any of claims 1 to 10. 

12. A coxoputer readable medium caanrying a contputer 

30 program product conprising elements adapted to carry out 
the method steps of any of claims 1 to 10 . 

13. Apparatus for analysing a heart beat signal which 
comprises a series of heart beat coxnplexeSe comprising t 

SS an extract set element (44) adapted to select a set 

of complexes from said series; 
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a select subset element (46, 48, 50) adapted to 
identify from the set a subset of complexes each, of which 
satisfy a similarity criterion with respect o members of 
the set; and 

a coiBbiner element (52) adapted to calculate a 
representative complex from the members of the subset. 

14. The apparatus of claim 13 further ccniprisiiig an 
analysis element (54) adapted to establish one or more 
heart beat signal parameters from said representative 
comolex. 
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15. The apparatus of claim 13 or 14 wherein the select 
subset element conprises: 

an align set element (46) adapted to align the 
con^plexes of said set; and 

a calculate criterion element (48) adapted to 
establish said similarity criterion from the aligned 

16. The apparatus of claim 15 wherein the calculate 
criterion element (48) establishes said similarity 
criterion as a threshold calculated using the conplexes of 
said set and an average of the members of the set. 

17. The apparatus of any of claims 13 to 16 further 
comprising a preprocess element (42) adapted to identify 
end select genuine heart beat c<mpl&z&B to form said 



18. Apparatus for analysing a heart beat signal 
comprising a plurality of heart beat con^lexes, the 
apparatus comprising: 

an analysis engine adapted to automatically calculate 
a representative heart beat cojnplex for each of a 
plurality of preselected intervals of said signal; and 

an editor adapted to enable a user to edit parameters 
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relating to each automatically calculated representative 
heart beat complex. 

19. The apparatus of claim 18 wherein each interval 
contains a set of heart beat con5>lexes, the analysis 
engine is adapted to automatically decide which complexes 
in the set to include in each representative heart beat 
GompleXi and the editor is adapted to enable the user to 
subsequently change which complexes axB included. 



20. 'The apparatus of claim 19 wherein the editor is 
adapted to enable the user to select whether a particular 
congplex should be included, not included, or included 
subject to a similarity criterion, in the corresponding 

15 representative heart beat congplex. 

21. The apparatus of any of claims 18 to 20 wherein the 
representative complex is an average of the included 
cos%)lexes . 

20 

22. The apparatus of any of claims 18 to 21 wherein the 
editor is further adapted to provide the user with a 
graphical display of at least some of the heart beat 
congplexes within a c\xrrent one of the intervals of said 

25 signal and with a graphical display of the current 
representative heart beat conqplex for the interval - 

23. The apparatus of claim 22 wherein the editor ie 
further adapted to display, on the display of the current 

30 representative heart beat complex, one or more feature 

markers indicating the automatically determined locations 
of one or more features of said representative heart beat 
complex. 

35 24. The apparatus of claim 23 wherein the editor is 

further adapted to enable the user to move said markers. 



« 



25. l»he apparatus of any of claims 18 to 24 wherein the 

analysis engine is further adapted to calculate one or 

more heart beat parameters of each representative heart 

beat coittplex. and said editor is adapted to display said 

5 heart beat parameters and to update said display according 
to changes made by said user. 

26. A con^uter readable medium con©rising a coxnputer 
program product arranged to provide the apparatus of any 

AO of claims 18 to 25 when executed on a suitable computer 
system « 
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There are disclosed methods and apparatus for 
aaalysing heart beat signals which comprise a series of 
similar heart beat coxnplexesj such as £CG signals. A 
subset of complexes satisfying a mutual similarity 

10 criterion is selected from a larger set. A representative 
complex is calculated from the members of the subset. One 
or more heart beat signal parameters are then established 
from the representative conplex. Computer apparatus 
making use of this scheme, including a graphical user 

15 interface for validating the representative complexes^ is 
al$o disclosed. 
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